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SUMMARY 


A parameter estimation technique is presented to estimate the radiative 
flux distribution over the Earth from radiometer measurements at satellite alti- 
tude. The technique analyzes measurements from a wide field of view (WFOV) , 
horizon to horizon, nadir pointing sensor with a mathematical technique to 
derive the radiative flux estimates at the top of the atmosphere for resolution 
elements smaller than the sensor field of view. A computer simulation of the 
data analysis technique is presented for both Earth emitted and reflected radia- 
tion. Zonal resolutions are considered as well as the global integration of 
plane flux. An estimate of the equator-to-pole gradient is obtained from the 
zonal estimates. Sensitivity studies of the derived flux distribution to direc- 
tional model errors are also presented. In addition to the WFOV results, medium 
field of view results are presented. 


INTRODUCTION 

The ability to predict climate is becoming increasingly important to the 
well-being of the United States and the world. The importance of the Earth's 
radiation budget in determining climate has long been recognized. The Earth- 
orbiting satellite provides a platform, outside the Earth's atmosphere, which 
is capable of simultaneously monitoring the outgoing reflection of the Sun's 
energy from the Earth's surface and atmosphere and monitoring the long-wave 
radiation emitted by the Earth and its atmosphere. Equally important, satel- 
lites are capable of monitoring these fluxes on a daily, monthly, or even a 
yearly basis. These capabilities provide the opportunity to conduct detailed 
studies of the variations in the radiation budget of the Earth, the effects of 
natural and man-made changes in the environment, and the effects which changes 
in the energy budget produce on the weather and climate of the Earth. 

The simplicity of the wide field of view (WFOV) radiometer makes it desir- 
able for Earth radiation mapping purposes if adequate spatial resolution can be 
obtained. The problem is depicted in figures 1 and 2. The WFOV radiometer 
receives radiation from all points on the Earth's surface from nadir to horizon. 
(See fig. 1.) For an orbit with an altitude of 600 km, the field of view is a 
circle with a diameter of 48° great circle arc. The radiation leaving the atmo- 
sphere is directionally dependent and thus so is the radiation impinging on the 
sensor at satellite altitude. Furthermore, as seen in figure 2, this direc- 
tional dependence is different at the two altitudes. The resolution which can 
be obtained with a set of overlapping measurements is better than that of a sin- 
gle measurement, as is illustrated in figure 3. For example, the location of a 
hot spot in an otherwise constant field can easily be determined to be within 
the shaded area. Recent theoretical developments (refs. 1 to 4) indicate that 
the spatial resolution associated with a set of WFOV measurements can be better 
than the instrument field of view, although the resolution obtainable does not 
match that obtainable by a narrow field of view scanning radiometer. Refer- 



ence 4 presents a statistical filtering technique for parameter estimation which 
inverts WFOV satellite data to obtain estimates of the radiative fluxes at the 
top of the atmosphere. The technique developed requires an a priori knowledge 
of the directional characteristics of exiting radiation. Because a satellite 
measurement of flux is necessarily indirect, directional modeling is required 
of any technique for obtaining flux values on any scale smaller than the globe. 
(For example, see ref. 5. ) 

In this paper a data analysis technique is simulated to demonstrate its 
feasibility for both Earth-emitted and Earth-reflected radiation. For the pur- 
pose of this paper, Earth-emitted radiation is defined as that part of the spec- 
trum between 5.0 and 50 ym and Earth-reflected radiation is that part of the 
spectrum between 0.2 and 5.0 ym. No assumption is made about the design of the 
sensor except the field of view, a flat-plate cosine response, and the ability 
to distinguish between the two wavelength bands. Mathematical models are devel- 
oped to represent simulated measurement data from a WFOV radiometer onboard an 
orbiting satellite. The simulated measurements are then analyzed by the param- 
eter estimation technique to yield an estimated radiation field at the top of 
the atmosphere. By comparing the estimated field to the "true" radiation field, 
the effects of directional model errors are studied. 

The parameter estimation technique is also simulated for a medium field of 
view (MFOV) radiometer. Here the field of view does not extend from nadir to 
the horizon but is restricted to a smaller area. Measurements from this sensor 
are simulated and then analyzed to yield an estimated radiation field at the 
top of the atmosphere. 


SYMBOLS 

A area of surface element, m 2 

a albedo 

B matrix of influence coefficients or bidirectional function 

ECA Earth central angle, deg 

F directional function for emitted radiation, sr~l 

f true anomaly of satellite, deg 

G equator-to-pole gradient, W/m 2 

g(a) angular response of sensor 

H solar constant (1350 W/m 2 ) 

h satellite altitude, km 

I radiation intensity, W/m 2 / /sr 
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inclination of orbit plane, deg 
number of measurements 

measured flux at satellite altitude, W/m 2 

square root of number of measurements, N 2 = J 

radiative flux at top of atmosphere, W/m 2 

mean global radiative flux at top of atmosphere, W/m 2 

radius of Earth-atmospheric system, km 

distance from surface element to satellite, km 

measurement system error, W/m 2 

cone angle from satellite nadir to point on surface of Earth, deg 

clock angle from North about satellite nadir to point on surface of 
Earth, deg 

incremental change 

Earth central angle between surface element and satellite, deg 

a specific value of Yr deg 

azimuth angle of exiting ray, deg 

variable of integration 

colatitude, deg 

directional model parameter 

zenith angle of Sun, deg 

bidirectional function, sr -1 

standard deviation, W/m 2 

variance of measurement system errors, W 2 /m^ 
longitude, deg 

zenith angle of exiting ray, deg 
longitude of ascending node, deg 
argument of periapsis or solid angle, deg or sr 


3 


Superscripts: 
m model 

T transpose 

t true 

Subscripts: 
max maximum 

o initial 

s satellite 

© Sun 

Nomenclature: 

FOVj field of view for jth zone 

WFOV wide field of view 

MFOV medium field of view 

A prime denotes a zonal value. A bar over a symbol denotes a vector and 
a circumflex C) over a symbol denotes an estimate. 


FORMULATION OF PROBLEM FOR EMITTED RADIATION 


The Earth atmosphere system will be approximated by a sphere of radius R. 
The intensity I of radiation leaving any point on this spherical surface is a 
function of colatitude 0, longitude <j>, the zenith angle 4 1 / and azimuth C of 
the exiting ray. (See fig. 4.) The emitted radiative flux Q at a point on 
the surface is given by (ref. 6) 


n TT/ 

=f 


I(0,<j>,\|;,5) cos sin ^ d^ dC 


A radiation directional function for emitted radiation F(0,cJ>,\p,^) called a 
limb-darkening function may be defined so that 


1 

i(0,4>,4»,£) = - Q(0,4>) F(0,<t>,\p,C) (2) 

7T 


In order that equations (1) and (2) be compatible, it is necessary that F 
satisfy the normalizing condition 
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F (0, C) cos sin \p d£ = 1 


(3) 


1 

7T 


n TT/2 
.'=0 


For this study F is assumed to have no azimuthal dependence- If the surface 
is Lambertian, that is, if radiation leaves the surface with equal intensity 
in all directions, then F(Q,<|>,^) = 1. 


Now consider a radiation sensor at satellite altitude, as shown in fig- 
ure 4. The radiation from the Earth incident on the sensor at a colatitude 
Q s , longitude 4> s , and altitude h is 

m(9 s ,<l> s ,h) = / i(0,4>,i|\S) g(0t) d<*> (4) 

Jfov 

where ^ and C are functions of the satellite position (0 S ,4> S ) and the sur- 
face element at the "top" of the atmosphere, w is the solid angle at the 
satellite subtended by the surface element, ot is the nadir angle at the 
satellite from the local vertical to the surface element, and the integration 
is carried out over the field of view FOV. The function g (ot) is the angu- 
lar response of the sensor to incoming radiation. For a perfectly black flat- 
plate sensor normal to the vertical (nadir looking), g(0t) = cos ot. Surface 
properties for other types of sensors can be incorporated in g(ot). Although 
the angular response function may depend on other variables such as azimuthal 
angle, this paper considers just the class of functions that depend on the nadir 
angle ot. By use of equation (2), equation (4) may be written as 


m(9 s ,<|> s ,h) 


1 

n 



Q(9,4>) F(0,<t>,^) g(a) du 


(5) 


Rather than integrate over the solid angle, it is convenient to integrate over 
the surface element dA or 


I f cos 

m(0 s ,4> s ,h) =- / Q(9,4>) F(9,<j>,40 g(a) dA (6) 

II Jfov r2 

where r is the distance from the satellite to the surface element. The objec- 
tive of the data analysis is to solve equation (6) for Q (9 r 4>) • 


The data analysis model of the radiation field will consist of 18 latitude 
zones, each 10° wide, in which the emitted flux is assumed to be constant. 

Thus, the measurements are modeled as 


1 18 

m (9-, cj> s , h) = — X Q-! 

IT j=l J 



F(0,<t>,Tp) g (a) 


cos t|> 
r2 


dA + v 


where Q-j is the emitted flux in the jth zone and the integration is over 
the surface area that is in both the field of view and in the jth zone. The 
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random error in the measurement is denoted by v. Thus, the ith measurement 
is modeled as 


18 

n»i - 23 Bi 3 Q j + Vi 

j=l 

where Bij is the influence coefficient of Qj on the ith measurement and 
is given by 


B 


ID 



F(0,4>,i|>) g (a) 


cos ^ 
r2 


dA 


(8) 


The variables a, and r are dependent on the satellite position, 0 S and 

<j) s , at the ith measurement. In matrix form, equation (7) becomes 


m = BQ' + v 


where m is a column vector of J measurements, Q' is a column vector of 
the 18 zonal fluxes, v is a column vector of J measurement errors, and B 
is the observation matrix with elements Bij. 

For a large measurement set (J » 18) , it is desired to find the best esti- 
mate of Q'. By the Gauss-Markoff theorem, the best (minimum- variance) linear 
unbiased estimate of Q' is given by the least-squares estimate (ref. 7, 
p. 349), that is, 

Q' = (B T B)" 1 B T m (9) 


where the measurement errors Vi are uncorrelated random variables with mean 
zero and variance CT^. Equation (9) together with the necessary supplementary 
relationship form the basis of the data analysis technique for estimating the 
emitted flux at the top of the atmosphere. 


SIMULATION PROCEDURE FOR EMITTED RADIATION 

The data analysis simulation consists of two basic parts, a true Earth- 
atmosphere radiation system for producing simulated data and a "model" Earth- 
atmosphere radiation system for the data analysis. A known radiation field is 
used to generate satellite measurements which are corrupted with random errors 
and passed to the model system as simulated data. The data analysis technique 
based on equation (9) is applied to these data and the zonal fluxes are esti- 
mated. The comparison is then made between the true flux and the estimated 
flux. This comparison indicates the accuracy of the data analysis technique. 

The simulation procedure as outlined above is graphically presented in figure 5. 
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Simulation of Emitted Radiation Data 


There are a number of concerns involved in generating simulated measure- 
ments. First, a typical radiation field must be defined over a grid system. 

The directional dependence F of the radiation must also be defined. In addi- 
tion, an orbit must be defined along with the number of measurements desired 
and the sampling strategy. Finally, a numerical integration scheme to obtain 
the measurements must be established. 

The grid system used in this study is illustrated in figure 5 and is an 
"igloo" grid system. Basically, the surface of the Earth is divided into 
squares of approximately equal area, the sides of which are latitude and merid- 
ian lines. For a 5° by 5° grid, the surface is divided into 5° latitudinal 
bands and each band is divided into an integral number of areas so that each 
area is approximately square. The band from 0° to 5° colatitude contains three 
areas. The band from 5° to 10° colatitude contains nine areas, and so forth. 
The total number of areas for a 5° by 5° grid is 1654. This is the grid system 
used throughout the study to simulate measurements. 

At each grid area the true emitted flux 4>i and its true directional 
dependence Fjr(4>) must be defined. The true emitted radiation field for 

this study was obtained from reference 5 which used the Nimbus III data. A 
map of the emitted flux is presented in figure 6. The directional functions 
used are the simple Lambertian function 

Lambertian: 



I 

fi 


F(40 = < 

1° 

and the 

empirical 

functions 

ence 5 , 

which are; 


Nominal 

: 

r 



1.074 exp 


F(^) = \ 

1. 074 exp 


0 

v. 


(0O < S go©) 

{y\) > 90°) 


(10a) 


(0° ^ < 60°) 
(60° ^ ^ S 90°) 
> 90°) 


(10b) 
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Upper 


1.03175 

1.03175 exp(o.l625{'-l + exp [(sec 30° - 
F(ip) = < ' 


(0° 2 i|> < 30°) 
sec iJ0/2.4lf}j 

(30° S \|i 2 90°) 
(\|» > 90°) 


(10c) 


Lower : 


F(\p) 


- 


1.106 exp^O. 0547 (1 - sec \|j) 


1 ° 


+ 0. 09375^exp[2. 24 (1 - sec «] - 1» 

(0° S \|i S 90°) 
(\p > 90°) 


(10d) 


where the functions are normalized by equation (3). The nominal function is 
a curve fit to a family of limb darkening functions based on Nimbus III data. 
The upper and lower functions represent upper and lower bounds, respectively, 
of the family of functions. 


The sampling strategy defines the location of the satellite for each mea- 
surement. The satellite orbit can be either elliptical or circular. Most of 
the numerical results are for a circular orbit with an altitude of 600 km. Fig- 
ure 7 defines the satellite geometry. The orbital inclination (i) was taken 
as 100° which corresponds to a common Sun synchronous orbit. The initial argu- 
ment of periapsis 0 ) o and the initial longitude of the ascending node ft 0 were 
taken as 180°. For this study periapsis of a circular orbit is defined arbi- 
trarily by the value given to 0) o . To state the sampling strategy simply, the 
measurements are uniformly spaced in true anomaly f and longitude. Most satel- 
lite measurements are equally spaced in time. The assumption of equal true 
anomaly spacing is equivalent to equal time spacing for circular orbits. This 
equivalence also holds approximately for orbits of small eccentricity. The uni- 
form spacing in longitude means that the ascending node is uniformly distributed 
between 0° and 360° after a sufficient number of orbits. The period of a typi- 
cal orbit in this study is about 1.5 hours so the satellite passes through 16 
ascending nodes a day. It will pass through 16 different nodes the following 
day, and so forth. Over the period of many days, the distribution of ascending 
nodes is well represented by the assumption of uniformity. 


Although in an actual mission there will be a large number of measurements, 
it is not necessary to simulate the total number. What is desired is to choose 
a representative set of simulated data which will allow the data analysis pro- 
cess to be studied. It was decided to generate measurements, where N 

measurements are computed on each of N revolutions of the satellite. Thus, 
the increment between data points is Af = 360°/N, and the increment in the 
orbit node after a complete revolution is Aft = 360°/N. The argument of peri- 
apsis is also incremented after a revolution by Aa) = Af/N so that measurements 
will be computed at a different set of colatitudes each orbit. The first revo- 
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lution is described by w = o) Q + Au) and ^ + Aft, and the kth revolution 

by 0) = o) Q + k Aa) and ft = ft Q + k Aft. in this manner, measurements are com- 
puted which are uniformly spaced in true anomaly and longitude. The measure- 
ment is then corrupted with a random error to simulate measurement errors. The 
errors are assumed to be normally distributed with mean zero and variance Cf 2 
and uncorrelated from measurement to measurement. 

The numerical integration scheme used to obtain the simulated data follows 
directly from the grid system. Since the true emitted flux field is defined by 
a 5° by 5° grid, the measurement (eq. ( 6 )) is approximated by 

1 ipt QiH Fi(^i) g(«i) cos i|»i 

m(9 s ,4> s ,h) = - > (11) 

i = l r i 

where the subscript i denotes each grid area. The zenith angle nadir 

angle o t^, and distance r:± between the satellite and grid area are computed 
for the center of the grid area defined by the mid-latitude and mid-longitude 
point. The field of view is defined by = 90° since = 90° is the hori- 
zon. If ipi < 90°, then the entire grid area is considered in the field of view 
and contributes to the measurement. If > 90°, then Ft-(^i) = 0, and this 

condition eliminates the grid areas not in the field of view. The measurement 
is then corrupted by adding a normally distributed random error. 


Data Analysis Model 


The true radiation system and the model radiation system are quite differ- 
ent as seen in figure 5. The modeled radiation field shown consists of 10° 
latitude bands each of which represents an area of constant flux. The direc- 
tional function is modeled as F m (^) and can be any of the directional func- 
tions given in equation (10). The superscripts on F m (40 and F t (4 J ) denote 
that the directional function may be different in the true system and the model 
system. Thus, the data analysis technique is applied to the measurements to 
estimate the 18 zonal fluxes which define the model system. The least-squares 
estimate of these fluxes is obtained from equation (9) where the elements of 
the B matrix are defined by the integral equation ( 8 ). The numerical integra- 
tion scheme used to evaluate these elements follows from the choice of the 
5° by 5° grid system? that is, the integral is approximated by 



Ai Fi(^i) g ( 04 ) cos 


r 


2 

i 


(12) 


where i ranges over all the grid areas that are in both the field of view 
and the jth latitude zone. The index k denotes the kth measurement and 
defines the position of the satellite 0 S and <|> s . Notice that the influence 
coefficient B^j is a function of the modeled direction function F m (^) 
whereas the simulated measurement (eq. ( 11 ) ) is a function of the true direc- 
tional function F t (4>). If F m (^) and F t (4 ; ) are different, then the esti- 
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mated flux will be in error and will indicate the sensitivity of estimates to 
the directional function. 


Properties of the emitted flux field which are important data analysis 
results are (1) the 10° zonal fluxes, (2) the global mean flux, and (3) the 
equator- to-pole gradient. Sets of simulated measurements were used to estimate 
these properties by using a variety of assumed directional models, F m OjO and 
FtftlOr and the results were compared with the true properties. In some cases 
F m (\p ) and Ft(ijJ) were the same and in other cases they were different. 

The estimate of the flux in the jth latitude zone is defined by Qj and 
is the jth element in the vector given by equation (9). This quantity is 

compared with the average true flux which is given by 


^ ] A iQi y y A iQi 



i 


where i ranges over all grid areas in the jth latitude zone. The area and 
flux of the ith grid area are denoted by A^ and Q*, respectively, and Aj 
is the area of the jth zone. The error in the jth zonal estimate is given by 
Qj - Qj and the average zonal error is defined by 


Average zonal error - 


18 


E abs loj - Qj I 

j=l 



(14) 


or in terms of percent the average zonal error is 


Average zonal error 



abs 




(15) 


This is an area weighted average since the zones have different areas. An equa- 
torial zone is much larger than a polar zone and the errors are weighted accord- 
ingly. Also, the absolute value of the error is averaged to eliminate cancella- 
tion. The average zonal error is a figure of merit for evaluating the data 
analysis technique. 


The global mean flux is another property of interest and is estimated by 
averaging the zonal estimates; that is. 
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( 16 ) 


Qg 


18 

2 

j=l 

18 

2*i 

j=l 


This estimate is compared with the true global mean flux given by 
18 1654 

S A i Q i ^ Ai ® i 

. = 

9 18 1654 

2 *5 2 m 

j=l i=l 


( 17 ) 


The definition of the equator-to-pole gradient adopted for this study is 
the average flux at the equation minus the average flux at the pole. Actually, 
the surface of the Earth is divided into four zones of equal area as shown in 
sketch (a) and the average flux for each zone defined. A northern and a 



60 ° 

= 90 ° 

= 120 ° 


southern gradient is defined and the two averaged to yield the equator-to-pole 
gradient. The average flux in the northern cap is an area-weighted average of 
six 10° zonal fluxes and is given by 
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On. cap 


(18) 


2 A j Q j 

j=i 


j=i 


The average flux for the other three zones is defined in the same manner and 
the gradient is 

(QN.eq ~ QN.cap) + (Qs.eq “ Qs.cap) 

G = 

2 


(19) 


The estimate of the equator-to-pole gradient is determined by replacing Q-j 
with Qj. 


RESULTS AND DISCUSSION FOR EMITTED RADIATION 

It is not practical to simulate the large amount of data of an actual mis- 
sion. The effect of large data sets is to reduce the error in the estimator as 
a result of improved spatial coverage and to reduce the variance in the estimate 
due to the measurement errors. This effect is illustrated in figure 8 where the 
average zonal error is plotted for various data sets. All measurements were 
made at a satellite altitude of 800 km and the true and model directional func- 
tions were Lambertian for all grid areas. As the number of measurements 
increases, the average zonal error decreases and asymptotically approaches the 
bias error in the estimation technique. For 3600 measurements, the average 
zonal error is approximately 0.5 percent. Although it is not necessary to sim- 
ulate large data sets, the number of measurements should be large enough to 
eliminate the effect of the location of the measurements. The pattern of the 
measurements is dictated by the uniform sampling strategy, but the location of 
the individual measurements depend on the node points of the various orbits 
defined by ftp. For example, if J = 625 (25 measurements per orbit times 

25 orbits), then the nodes are approximately 14° apart. Thus, = 180°, 

183.5°, 187°, 190.5° implies four different measurement sets each containing 
625 uniformly spaced measurements. The variation in the average zonal error 
for these four cases is presented in figure 8. Since this variation due to 
sampling is unacceptably large, 1296 (36 2 ) measurements were chosen as an appro- 
priate data set for the simulation results presented herein. The reduction 
in sampling variation achieved by this is also shown in figure 8. 

The random errors added to the computed measurements were all from a nor- 
mal distribution with mean zero (unbiased) and variance 1 (a = 1 W/m 2 ) . This 
corresponds to a three standard deviation error of 3 W/m 2 or a 1-percent maxi- 
mum error when the measurement is 300 W/m 2 . This is a realistic accuracy for 
present day WFOV radiometers (ref. 8) . The errors in the flux estimates due 
to unmodeled longitudinal variation of the radiation field, however, far exceed 
the estimation errors due to random measurement errors. To see this effect, 
consider the computed measurements on an equatorial orbit. If the zonal flux 
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were constant , the standard deviation in the measurement would be near 1 W/m 2 . 
The actual standard deviation, however, is 19.7 W/m 2 as a result of the longitu- 
dinal variation of the true radiation field. Thus, the relative effect of ran- 
dom measurement errors or the simulation results was small. 

One of the major advantages of a simulation study, as compared with an 
error analysis, is that modeling errors can be investigated more easily. The 
model radiation system consists of 18 latitude zones of constant flux, which 
is very different from the true radiation system which consists of 1654 areas 
of constant flux. This spatial modeling error introduces bias errors into the 
estimation process. Another modeling error is encountered when F^(^) and 
F m (^) are different. The data analysis simulation demonstrates the accuracy 
of the data analysis in the presence of these modeling errors. 

Three simulations were performed to determine the estimation errors result- 
ing from the 10° zonal spatial model. In each case, the true and model direc- 
tional functions were the same for all grid areas, that is, Ft(\p) = F m (\p). The 
estimation errors, therefore, resulted only from the spatial modeling error and 
the random measurement errors. The orbit assumed for these three cases was cir- 
cular with h = 600 km, i = 100°, and 0) o = = 180°. The sampling strategy 

assumed was uniform consisting of 1296 measurements based on a flat-plate sen- 
sor (g(ot) = cos a) . in the first case F fc (^) and F m (\p) were both Lambertian 
(eq. (10a)). The second case assumed Ft(^) and F Itl (4 J ) were the nominal func- 
tions (eq. (10b)). In the third case, part of the surface radiated according 
to the nominal function and the remainder of the surface was Lambertian. The 
results from all three simulations were similar. 

The second simulation corresponding to the nominal functions is presented 
in figure 9. The 10° latitude zone estimates of the flux are plotted as a func- 
tion of latitude and compared with the true zonal flux (eq. (13)). These esti- 
mates represent the actual distribution very well. Furthermore, when used to 
estimate the global mean (eq. (16)), the error is less than 0.01 percent. This 
increase in accuracy results because the resolution is better for large areas 
than for small areas and because a flat-plate WFOV sensor is well suited to mea- 
sure the global average. In summing to get a global mean, many of the errors 
in the zonal estimates cancel as can be seen by the alternating signs on the 
errors. The underlying reasons for these alternating signs and the cancella- 
tion are discussed in reference 4. The high degree of accuracy for the global 
mean results partially because the error in the measurements is unbiased. If 
the measurements were biased, then the global mean would be biased accordingly. 
The error in the equator-to-pole gradient was 1.22 percent. 

These three cases demonstrated that the spatial resolution obtained from 
a set of WFOV measurements is smaller than the instrument field of view and 
that the flat-plate WFOV sensor is well suited to measure both zonal fluxes 
as well as the global mean flux. The zonal flux estimates, based on zones of 
constant flux, represent the average of the true fluxes over a zone. This 
results from the nature of the WFOV sensor and the uniform sampling strategy. 

The estimation process worked equally well for uniform directional functions 
(Lambertian) and nonuniform directional functions (limb darkening) . Defining 
different directional functions for different areas of the surface caused no 
problem in estimation. 
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A fourth simulation similar to the first (Lambertian) case was performed 
to test the sensitivity of results to orbit parameters. This case used an ellip- 
tical orbit with a periapsis altitude of 600 km and an apoapsis altitude of 
800 km. Accuracies obtained were similar to those of the first three cases. 

The effect of directional function model errors on the emitted flux esti- 
mates was also studied. When F^M) and F m (xp ) are different, the result is 
a bias of the zonal estimates. These biases do not become negligible with an 
increasing number of measurements. Several simulations were performed to study 
these biases for different directional functions and the results are presented 
in table I. In the first simulation, the true directional function is assumed 
to be the nominal function and the model function is Lambertian. This simula- 
tion represents the situation where the directional dependence of the radiation 
field is ignored and is assumed to be uniform in the data analysis. For this 
case with 1296 measurements the average zonal error is biased by 0.70 percent 
or 1.67 W/m^ and the global mean error remains small. The equator-to-pole gra- 
dient is biased 1.65 percent. The next two simulations represent a more rea- 
sonable approach to modeling which assumed the nominal function for F m (^) 
whereas F^- (4*) was taken as the upper and lower limb darkening function, respec- 
tively. If one assumes that these simulations establish reasonable bounds on 
the bias error due to directional function modeling, then these two cases indi- 
cate that the average zonal error will be less than 1.37 W/m^ if the radiation 
is modeled by the nominal limb darkening function. The next simulation was pre- 
sented previously in figure 9 and showed that perfect directional modeling 
(F fc ) = F m (40) results in a 1.30-W/m^ zonal error. This error is mainly due 
to the spatial modeling error which introduces a bias in the estimates; that is, 
the longitudinal variation in the radiation field is not modeled. The random 
error in the measurements also contributes to this error. For example, if the 
true radiation field were constant within a zone, then the average zonal error 
would be reduced from 1.30 W/m^ to 0.24 W/m^. This remaining error is the 
result of measurement error alone and approaches zero with increasing measure- 
ments at the rate l/\/j. In addition, if the measurements were error free, 
then the average zonal error would be zero. The worst case which compares the 
true lower function and the model Lambertian function has an error of 1.76 W/m^ 
which is only 0.46 W/m^ above that which could be expected from perfect direc- 
tional modeling. Finally, table I shows that an elliptical orbit introduces 
no additional problems. 

Three error sources have been defined and studied. The spatial modeling 
error and the directional modeling error cause biases in the zonal flux esti- 
mates. The random measurement error causes a random error in the zonal flux 
estimate which tends to zero with increasing number of measurements. The errors 
caused by these three error sources interact with each other to produce the 
final results. This interaction must be considered when interpreting the 
results of table I. Nevertheless, these simulations imply that the error asso- 
ciated with modeling the radiation field as 18 latitude zones of constant mag- 
nitude flux is greater than the error associated with modeling the directional 
dependence. This implies that a finer resolution spatial model would be bene- 
ficial for obtaining zonal flux estimates. 
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II 


FORMULATION OF PROBLEM FOR REFLECTED RADIATION 


The major difference between the formulation of simulations for emitted 
radiation and reflected radiation is the inclusion of the Sun as the radiation 
source for reflected radiation. The intensity I of radiation reflected by 
any point at the top of the atmosphere is a function of colatitude 0 , longitude 
<t>, the zenith angle 4* of the reflected ray, the zenith angle 5 of the Sun, 
and the azimuth £ of the reflected ray relative to the Sun. The reflected 
flux Q at the point is given by 

J f2fr rr 1/2 

I cos 4 > sin 4 > d £ ( 20 ) 

z=o Ji>=o 

A radiation directional function p ( 0, 4>, \p, £,£) , called a bidirectional reflec- 
tance function, may be defined by 

I(0,4>,4>,C,S) = Q(9»4>f ?) P(0,4>,^f CrO (21) 


In order that equations (20) and (21) be compatible, it is necessary that p 
satisfy the normalizing condition 


n TT/2 
.'=0 


P(9,<t>,4v C,£) 


cos 4 1 sin 4 1 d 4 > d£ = 1 


( 22 ) 


Also, the reflected flux can be modeled as 

f H a( 0 , 4 >) cos £ 

Q (0,<t>,Z) = 


(0° ^ ^ 90°) 
(90° < £ ^ 180°) 


(23) 


where H is the solar constant and a (B r cj>) is albedo at the top of the atmo- 

sphere. Note that albedo is assumed to vary spatially but not temporally. Thus, 
the diurnal effect of reflected flux is modeled by cos £ and allows for varia- 
tion only due to the rising and setting of the Sun. The radiation data then 
will be analyzed by the parameter estimation technique to produce an albedo map 
at the top of the atmosphere which is related to reflected flux by equation (23) 
depending on the Sun's location. To proceed, the measurement is given by 


m ( B 3 , <J>g, 0Q, 4*0/ h) 



I(9/<t>f^fC,5) g(ot) do> 


where a) is the solid angle at the satellite subtended by the surface element 
at the top of the atmosphere and 0 q and 4 )q denote the position of the Sun. 
Substitution of equations (21) and (23) yields 


m( 0 s , (J> s , 0 0 , 4 > 0 , h) 



H a (0 , 4>) cos Z, P (0,<|>,iK £,5) g(a) dw 
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or 


m(0, 


5» <t > S'®0' < t ) 0'^) = / 

Jfov 


cos 4* 

H a(0,4>) cos 5 p(0,<MrC»5> g(«) dA 


where r is the distance from the satellite to the surface element dA. The 
data analysis model of the albedo field will consist of 18 latitude zones/ each 
with constant albedo aj. Thus, the measurements are modeled as 


m(0< 


, f cos 4) 

/<f>g, 00,4>0, h) = aj I H cos E, p (0 ,4* ,4 * fK ) g(ot) 

i=i Jfov* r2 


dA + v 


where the integration is over the surface area that is in both the field of view 
and the jth zone. The random error in the measurement is v. Thus, the ith mea- 
surement is modeled as 
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m 


i “ ^ B ij a j + v i 
j=l 

where B-h is the influence coefficient of a-j on the ith measurement and 
is given by 


B 


ID 


■/ 

JFOVa 


cos 

H cos £ p (Q'ty'tyiZ'K) 9(a) dA 


(24) 


In matrix notation 

m = B5' + v 


where m is a column vector of J measurements, 5' is a column vector of 
the 18 zonal albedos, v is a column vector of J measurement errors, and 
B is the observation matrix with elements B|j (eq. (24)). According to equa- 
tion (9) , the least-squares estimate of albedo is 


a’ = (B T B)" 1 B T in 


(25) 


and the resulting estimate of flux is given by equation (23) 


SIMULATION PROCEDURE FOR REFLECTED RADIATION 

The simulation procedure for reflected radiation is very similar to the 
procedure for emitted radiation. Consider a true Earth-atmosphere radiation 
system and a model Earth-atmosphere radiation system and compare the true flux 
with the estimated flux. The estimated parameter, however, is albedo which is 
related to flux through dependence upon the position of the Sun. 


16 



The true system consists of a 5° by 5° grid system of 1654 grid areas. 

The true albedo for this study was obtained from reference 5 which used the 
Nimbus III data. 

The bidirectional function Pf'KCr?) has been assumed to vary only with 
the zenith of the Sun f the zenith of the reflected ray, and the azimuth between 
the two and does not vary with latitude or longitude. A sketch (sketch (b) ) of 
the geometry at a point P at the top of the atmosphere shows the three vari- 
ables where 

OO < \\) ^ goo 

0° ^ K = 90° 

0° ^ e ^ 360° 

The notation "forward" denotes the forward reflections of the Sun. 


Sun Vertical 



Sketch (b) 

The functional form for ^ = 90° is given by 

= — - — /l + Ci (sin + C 3 sin 2 / 2 ^ + C 4 sin® £ sin 2 (cos 2 £ 


+ C 5 sin 2 S) 


Cg cos ?\/l - sin 2 C sin 2 £ 


s in 


4> 1 


a 


+ -(1 + 3 cos 4?) 
2 



(26) 
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where 


C! = 1.5853 
C 2 = 13.4646 
C 3 = -0.01717 
C 4 = -1.856 
C 5 = 0.3046 
C 6 = -0.357 

and N(£) is the normalizing function so equation (22) is satisfied. A graph 
of N(£) is presented in figure 10 and p(\|;,£,£) is presented in figure 11. 
The bidirectional function is set equal to zero for ^ > 90°. This bidirec- 
tional function is an empirical model based on the data of reference 5. To 
simulate bidirectional modeling error, a family of functions has been defined 
and is given by 


1 

Boi>,e,£) = x poi>,c,S) + (i - X)- (oa^ d ( 27 ) 

TT 

The bidirectional function B is a composite of the empirical model and a dif- 
fuse model and represents a linear variation with X between the two models. 
Thus, to simulate bidirectional modeling error, one can generate the measure- 
ments with X = 1 (empirical model) and estimate the albedo with X = 0 (dif- 
fuse model) or any other combination of X values. The error in the estimate 
of albedo will then indicate the relationship between modeling errors and albedo 
estimates. 

For all the simulations presented here, the Sun's position has been assumed 
to be fixed with respect to the orbital plane. If the orbit were Sun synchro- 
nous, then this would be the case and the change in the orbit node with respect 
to the Sun would be zero. For low inclination orbits, this assumption is valid 
only for a few days. Since the true albedo field is constant with time (that 
is, it has no diurnal variation) , the diurnal variation of the flux field is 
deterministic and a result of the Sun elevation alone. Thus, data at one Sun 
angle or one local time is sufficient to estimate the albedo field. Fixing the 
Sun position relative to the orbital plane, therefore, is consistent with the 
assumption of a non- time-dependent albedo. Nevertheless, different Sun posi- 
tions are considered to determine their effect upon the albedo estimates. 

The numerical integration scheme used to obtain the simulated reflected 
flux measurement is based on the 5° by 5° grid system and is given by 

1654 H cos p(\J)i,Ci,5i) 9(«i) cos \|>i ai 

m(9s' = I ~ (28) 

<_ i r2 
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where the subscript i denotes each grid area, the bidirectional function has 
been assumed equal to P for all grid areas, and ax is the true albedo for 
each grid area. The wide field of view is defined by £ 90°. The 5° by 5° 
grid system is also used to approximate the influence coefficients (eq. (24) ) ; 
that is. 


(29) 

* i 

where i ranges over all the grid areas that are in both the kth field of view 
and the jth latitude zone. The influence coefficients reflect directional model- 
ing errors which are controlled by the parameter X. 


B kj - J2 


H cos g(<*i) cos VjAi 


(1 - X) 

Xp(i|»i,Ci,Ci) - 

TT 


The estimates of the zonal albedos are given by equation (25) and they are 
compared with the true zonal albedos given by 

23 A i a i 


J L, A i 

i 

With the position of the Sun known, one can convert the average zonal albedo 
into the average zonal flux. This relationship is derived in the appendix. If 
the Sun is in the equatorial plane, then the average zonal flux is 



ajH(A0 - sin A0 + 2 sin A0 sin 2 0j) 


4tt sin 



sin 0j 


(31) 


where 5j is the mid-latitude of the jth zone or = 0j + (A0/2). The 

estimates of zonal flux are obtained by evaluating equation (31) with a j = a j- 
These fluxes are then used to evaluate the average zonal error (eqs. (14) and 
(15)), the global mean flux (eqs. (16) and (17)), and the equator-to-pole gradi- 
ent (eqs. (18) and (19)). 


RESULTS AND DISCUSSION FOR REFLECTED RADIATION 

The model radiation system for the reflected simulation study consists of 
18 latitude zones of constant albedo. The results from a typical simulation 
with no directional modeling errors are presented in figure 12. The estimate 
of the albedo field is seen to represent the true field very well. Also shown 
is the reflected flux for the Sun in the equatorial plane. The odd number of 
measurements results from considering 2500 uniformly spaced measurements, 917 
of which were in total darkness. Thus, only 1583 measurements were used to 
estimate the zonal albedos. The average zonal flux error for this simulation 
was 1.24 W/m 2 and the global mean error was -0.17 W/m 2 . The equator-to-pole 
gradient was in error by -0.49 W/m 2 . Since the gradient is a difference between 
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fluxes, its error is not reflected directly by the zonal or global errors. For 
example, if the albedo estimates were uniformly high, then the zonal and global 
errors would be affected directly while the gradient would only be indirectly 
affected. Suppose the situation was such that the zonal errors canceled each 
other and the global error was zero. If the zonal errors were evenly distrib- 
uted over the globe, then the gradient error would be minimal. However, if the 
zonal estimates were high in the polar regions and low in the equatorial region, 
then the gradient error would be much larger than one would expect by consider- 
ing only the zonal and global errors. 

The directional modeling error has been studied by assuming a bidirec- 
tional model that varies linearly with X between the empirical bidirectional 
model and a diffuse model. The simulated measurements were computed by using 
p(4>»£,£) (eq. (28)) and the albedos were estimated by using 

1 

B(\|»,c,5) = Xp(\|),C,5) + (1 - X)- 

7T 

as the directional model (eq. (29)). The diffuse assumption corresponds to 
X = 0 and perfect directional modeling corresponds to X = 1. Results of 
simulations for various X values are presented in figure 13. As would be 
expected, the global and zonal errors are reduced with better directional mod- 
eling. They exhibit errors, however, even for perfect directional modeling. 

These errors are the result of the spatial modeling error; that is, the error 
associated with representing the albedo field by 18 zones of constant albedo. 

The gradient error exhibits a behavior much different from the global or zonal 
errors. It has a positive error for X = 0 and a negative error for X = 1. 

For the reflected radiation simulations discussed, the position of the Sun 
was assumed to be at a right ascension and declination of zero; that is, the 
Sun was assumed to be in the orbital plane and at the equator. This assumption 
corresponds to a noon Sun synchronous orbit during equinox. The effect of chang- 
ing the right ascension of the Sun to 45° is presented in figure 14. This case 
is the same as the case presented in figure 12 except 4 >q = 45°. The individual 
zonal errors are plotted against latitude for Sun right ascension of zero and 
45°. The conclusion is that the right ascension of the Sun does not signifi- 
cantly change the results of the simulation. Hence, the right ascension of the 
Sun was assumed to be zero. If the declination of the Sun is changed, one of 
the poles will be in darkness and an albedo estimate of that pole is impossible. 
Six months later, however, the situation is reversed and the other pole is in 
darkness. Thus, over a 6-month period both poles are observable. For this 
reason the declination of the Sun was assumed to be zero throughout this study 
so that both poles would be observable. 


RESULTS AND DISCUSSION FOR MEDIUM FIELD OF VIEW SENSOR 

All the simulation results presented thus far have been for a wide field 
of view radiometer which receives radiation from all points on the Earth's sur- 
face from nadir to horizon. If the field of view is restricted, it becomes a 
medium field of view (MFOV) radiometer. The analysis of this instrument using 
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the simulation procedure outlined previously is straightforward. The only 
change is that a smaller surface area is viewed by the instrument. Since the 
sensor is still a nadir pointing instrument, its field of view can be defined 
by either a the cone angle from satellite nadir to the edge of the field of 
view or by Y the Earth central angle (fig. 4) . Previously, the field of view 
was defined by = 90°. The relationship between these two angles for various 
satellite altitudes is presented in figure 15. The curve labeled y = y* cor- 
responds to the WFOV sensor and at an altitude of 600 km, Y* = 23.93° and 
a = 66.06°. Of course, Y* changes for different altitudes. The other two 
curves represent a constant Y- For example, if Y = 5° and h = 600 km, then 
a = 41.68°. Another notation used to describe the MFOV sensor is the total 
Earth central angle (ECA) of the field of view or ECA = 2 y- All the simula- 
tions reported in this section are for an ECA of 10°. Results from a typical 
MFOV simulation for emitted radiation are shown in figure 16. Directional mod- 
eling errors were not considered for this case. It can be seen that the esti- 
mated flux field matches the real field well except at the North Pole. The 
average zonal error is nearly the same as for the perfect modeling WFOV case of 
figure 9. The global mean error and the gradient error, however, are larger for 
the MFOV case. When the directional modeling errors were considered, errors in 
the estimate of the zonal flux, global mean, and gradient increased even more 
significantly. The results are presented in table II. Recall that the WFOV 
flat-plate sensor is well suited to measure the global mean flux independent of 
the directional function. This is not the case for MFOV sensors. In fact, 
since the MFOV sensor does not receive radiation from high zenith angles, the 
global estimate is biased depending on the assumed directional functions. These 
biases also affect the accuracy of the zonal estimates. An interesting result 
is that the equator -to-pole gradient can be estimated better with a WFOV sensor 
than with a MFOV sensor. This condition possibly could result from the reduced 
spatial sampling of the MFOV sensor. 

The same type of analysis for the estimation of the albedo field was per- 
formed for the MFOV as for the WFOV sensor. The results of a typical simula- 
tion are presented in figure 17 and are very similar to the WFOV sensor results 
(fig. 12) . The sensitivity of these estimates to directional modeling errors 
for reflected radiation is shown in figure 18. These errors are generally 
larger than the WFOV errors except for perfect directional modeling (X = 1) . 
Thus, the WFOV and the MFOV estimates are comparable for perfect modeling. 
However, the MFOV sensor is much more sensitive to directional modeling errors 
than the WFOV sensor. 


CONCLUDING REMARKS 

A parameter estimation technique has been defined to analyze wide field of 
view (WFOV) radiometer measurements to obtain radiation flux estimates at the 
top of the atmosphere. A computer simulation of the data analysis technique 
has been performed for both Earth-emitted and reflected radiation. Results 
were also obtained for a medium field of view (MFOV) radiometer. 

The error sources considered in estimating the emitted flux at the top 
of the atmosphere for WFOV measurements are the assumed limb darkening func- 
tion, the assumed spatial model of the flux field, and the random error in the 
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measurement. If the flux field is modeled as 10° zonal bands of constant flux 
and a nominal limb darkening function is chosen, the average zonal error is 
1.37 W/m 2 . Most of this error is attributed to the spatial modeling. A per- 
fect limb darkening model only reduces the zonal error to 1.30 W/m 2 . If the 
zones are in fact regions of constant flux, the zonal error is reduced to the 
measurement error which is 0.24 W/m 2 for a sample size of 1296 measurements and 
decreases to zero with increased measurements. The spatial modeling error and 
the limb darkening error, however, result in an estimate bias that does not 
decrease to zero with increased measurements. The global mean error was found 
to be insensitive to the spatial and limb darkening models and was less than 
0.05 W/m 2 . The equator-to-pole gradient can be estimated to within 1.0 W/m 2 . 

The error sources considered for estimating the reflected flux at the top 
of the atmosphere from WFOV measurements are the same as those for the emitted 
flux. If the spatial model is assumed to be 10° zonal bands of constant albedo 
and the directional dependence of the flux is taken to be diffuse instead of 
bidirectional, then the average zonal error is 2.92 W/m 2 , the global mean error 
is -1.44 W/m 2 , and the gradient error is 0.50 W/m 2 . Removing the directional 
modeling error by choosing the correct bidirectional function reduces these 
errors. For perfect directional modeling, the zonal error is 1.24 w/m 2 , the 
global mean error is -0.17 W/m 2 , and the gradient error is -0.49 W/m 2 . 

Results have been obtained for a medium field of view radiometer. Instead 
of viewing all points on the surface from nadir to the horizon, the MFOV is 
restricted to an Earth central angle of 10°. For perfect directional modeling, 
the WFOV and MFOV estimates of the zonal and global flux are approximately the 
same. However, the flux estimates from MFOV measurements are more sensitive to 
directional modeling errors than are the flux estimates from WFOV measurements. 
For example, the emitted flux global mean could be in error by as much as 
5.55 W/m 2 for MFOV as compared with 0.05 W/m 2 for WFOV. This error also affects 
the zonal error for MFOV. The emitted flux equator-to-pole gradient can be 
estimated better with a WFOV sensor than with a MFOV sensor. This is in part 
a consequence of the reduced field of view and the large flux errors at the 
poles exhibited by the MFOV estimates. Results for the MFOV reflected flux 
estimates are also very sensitive to directional modeling errors. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 29, 1978 
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APPENDIX 


AVERAGE REFLECTED FLUX OVER AN AREA 
Consider the average reflected flux over an area A given by 


X Q dA I aH cos £ dA 

/A 


where a is albedo, H is the solar constant, and £ is the zenith of the 
Sun. For a latitude zone defined by 9i^0^0i+A0 and 0 ^ <}> < 2tt , write 


/ /*4>=2n /*9=0 i +A0 

dA= / I 

J >=0 ^ 0 =0 


R 2 sin 0 d0 dcf) 


A| = 4ttR 2 sin f— j sin 0^ 


where 0^=0^+ (A0/2) is the average colatitude and R is the radius of 
the Earth-atmospheric system. Before the average flux can be computed from equa- 
tion (Al) , one needs an expression for cos £ in terms of <j) and 0. The unit 
vector to the Sun is 

A 

Xq = (sin 00 cos <|> 0 , sin 0 @ sin 4 > @ , cos 0 @ ) 
and the unit vector normal to the surface element dA is 


X A = (sin 0 cos 4>, sin 0 sin cos 0) 


Thus, 


5 = *0 * x A 


= sin 0 q sin 0 cos <J>q cos 4 > + sin 0 ~ sin 0 sin sin + cos 0 ^ cos 0 


= sin 0 q sin 0 cos (<J) - 4>q) + cos 0^ cos 0 
and for a zone of constant albedo a', one has 


1 fa-H 
A ’ Ja 


cos £ dA 
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, /■<j>=2rr /*e=e+(Ae/2) 


(A0/2) 


H(sin sin 0 cos (4> 


+ cos 0 q cos 0)R 2 sin 0 <30 <34> 

Since the reflected flux is zero over the dark side of the Earth {£ > ir/2), 
one can redefine the variable of integration as D * <j> - <j>Q and integrate from 
-n* to n* where n* represents the Sun terminator which is considered con- 
stant over A0; that is. 


a . f n=n* T0=0+(A0/2) 

- 

-'T]=-D* ■'0=0- 


H(sin 0-. sin 0 cos ti + cos 0^ cos 0)R 2 sin 0 d0 dn 


(A0/2) 


Q' = — HR 2 (A0 - sin A0 + 2 sin A0 sin 2 0) sin n* sin 0, 
a» L '■ 


+ 2ri*(sin A0 sin 0 cos 0) cos 0 qJ 
where ti* is defined by cos £ = 0 or 

sin 0/-. sin 0 cos n* + cos 0^ cos 0 = 0 


n* = cos-1 (cot Q C ot 0) 


(0 ^ n* ^ TO 


If the Sun is in the equatorial plane, then = tt/ 2, n* = tt/ 2, and 


Q' = — HR 2 (A0 - sin A0 + 2 sin A9 sin 2 0) 

A' 


and from equation (A2) 


a'HR 2 (A0 - sin A0 + 2 sin A0 sin 2 0) 
4 ttR^ sin I — J sm 0 


a'H(A0 - sin A0 + 2 sin A© sin 2 0) 

/A9\ I 
4tt sin I — J sin 0 
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APPENDIX 


It is interesting to note that as the width of the latitude zone becomes small 
such that sin AQ ~ A0, the average zonal flux is 

a'H sin 0 Qmax 

q* * = 

7T 7T 

or the average zonal flux is approximately equal to the maximum flux at the 
Sun's meridian divided by tt which accounts for the reduction in flux away 
from the Sun's meridian. 
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TABLE I.- SENSITIVITY OF WFOV ESTIMATES TO DIRECTIONAL FUNCTION MODEL FOR EMITTED RADIATION 



orbit; 


1296 measurements; sensor error, 



F r 010 

F m (\p) 

Global mean error 

Average zonal error 

Gradient 

error 

Altitude, 

km 

Remarks 

Percent 

W/m 2 

Percent 

W/m 2 

Percent 

w/m 2 

Nominal 

Lambertian 

-0.0004 

-0.001 

0.70 

1.67 

1.65 

0.81 

600 

No model 

Lower 

Nominal 

-.0094 

-.022 

.57 

1.37 

1.16 

.57 

600 

Lower bound 

Upper 

Nominal 

.0188 

.045 

.45 

1.09 

1.07 

.52 

600 

Upper bound 

Nominal 

Nominal 

.0073 

.017 

.54 

1.30 

1.22 

.60 

600 

Perfect model 

Lower 

Lambertian 

-.0174 

-.042 

.74 

1.76 

1.58 

.78 

600 

Worst case 

Nominal 

Lambertian 

.0085 

.020 

.68 

1.62 

2.15 

1.06 

600 x 800 

Elliptical orbit 


TABLE II.- SENSITIVITY OF MFOV ESTIMATES TO DIRECTIONAL FUNCTION MODEL FOR EMITTED RADIATION 

[jL00 o orbit; 1296 measurements; sensor error, 1 W/m 2 J 


F r ty) 

F m (40 

Global mean error 

Average zonal error 

Gradient 

error 

Altitude, 

km 

Remarks 

Percent 

W/m 2 

Percent 

W/m 2 

Percent 

W/m 2 

Nominal 

Lambertian 

5.40 

12.94 

5.56 

13.31 

7.72 

3.79 

600 

No model 

Lower 

Nominal 

1.04 

2.49 

1.27 

3.05 

3.04 

1.49 

600 

Lower bound 

Upper 

Nominal 

-2.32 

-5.55 

2.41 

5.76 

-.11 

-.06 

600 

Upper bound 

Nominal 

Nominal 

.12 

.29 

.60 

1.43 

2.25 

1.11 

600 

Perfect model 

Lower 

Lambertian 

6.37 

15.26 

6.52 

15.62 

8.55 

4.19 

600 

Worst case 

Nominal 

Lambertian 

5.34 

12.79 

5.50 

13.17 

8.02 

3.93 

600 x 800 

Elliptical orbit 





Satellite 


Figure 1 



- Wide field-of-view radiometer data analysis problem. 




(a) Radiation at satellite altitude. 



(b) Radiation at top of atmosphere. 

Figure 2.- Radiation directionality at top of atmosphere and 

at satellite altitude. 



Orbit track 


Figure 3,- Information from measurements with overlapping view fields. 
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Figure 4,- Earth-satellite geometry. 
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Figure 5.- Simulation of emitted wide f ield-of-view 
measurements and analysis. 
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Figure 6.- 


u> 

u> 


Earth emitted flux map used in simulation 





Emitted zonal flux, W/m 


Sensor WFOV 
1296 measurements 
100° orbit - 600 km 
Uniform sampling 2 

Measurement noise 1 W/in 



90N 60 30 EQ 30 60 90S 


Latitude, deg 


Figure 9.- WFOV estimates of average zonal emitted flux. 
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Average reflected zonal flux, W/m 
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Figure 13.- Sensitivity of WFOV estimates to directional modeling errors. 
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i = 100°; h = 600 km; J = 1600; and a = 0. B = Xp - (1 - X)—. 
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Figure 14.- Effect of sun longitude on WFOV estimates of reflected flux, 
i = 100°; h = 600 km; X = 1; and sun latitude, 0°. 
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Figure 17.- MFOV estimates of average reflected zonal flux. 
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(No modeling) (Perfect modeling) 

Figure 18.- Sensitivity of MF0V estimates to directional modeling 
errors. i = 100°; h = 600 km; ECA = 10°; J = 1300; and a = 0. 
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B = Xp - (1 - X)-. 
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